Observation of Self-binding in Monolayer 3 He 
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We report clear experimental signatures of the theoretically unexpected gas-liquid transition in 
the first three monolayers of 3 He adsorbed on graphite. The transition is inferred from the lin- 
ear density dependence of the 7-coefficient of the heat capacity measured in the degenerate re- 
gion (2 < T <80 mK) below a critical liquid density (p c o)- Surprisingly, the measured p c o values 
(0.6~0.9 nm -2 ) are nearly the same for all these monolayers in spite of their quite different environ- 
ments. We conclude that the ground-state of 3 He in strict two dimensions is not a dilute quantum 
gas but a self-bound quantum liquid with the lowest density ever found. 

PACS numbers: 67.30.hr, 67.30.ej, 67.10.Db, 67.30.cf 



Matter can in principle be in either gas or a liquid 
phase at absolute zero if the quantum parameter, the 
zero-point kinetic energy divided by the potential en- 
ergy, is large enough. The two-dimensional (2D) helium-3 
( 3 Hc) system has long been thought as the only material 
which stays gaseous at the ground state [l[ . This system 
is experimentally realized in 3 He monolayers adsorbed on 
an atomically flat and strongly attractive graphite sur- 
face. Mostprevious theories based on the variational cal- 
culations [i-Q, the diffusion Monte Carlo calculation [H[ 
and the Fermi hypernetted chain method Q support the 
absence of self-binding of 3 He in 2D. Indeed, no signature 
of the gas-liquid (G-L) transition was experimentally ob- 
served in the first and second layer 3 He on graphite down 
to T ss 3 mK and to areal density p = 1 nm -2 0. This is 
in sharp contrast to monolayer 4 He with smaller quantum 
parameter on graphite. It is well established experimen- 
tally [|| and theoretically [i| that in this system the G-L 
transition takes place at temperatures below 1 K and the 
self-bound liquid density at T = (p c o) is 4 nm -2 . 

The first experimental address to this problem was 
made by Bhattacharyya and Gasparini [Io| , who found a 
kink or small discontinuity near 100 mK in the heat ca- 
pacity (C) of submonolayer 3 Hc floated on a thin super- 
fluid 4 He film adsorbed on a Nuclepore substrate. They 
attributed this to a puddle formation of 3 He in 2D. It 
is to be noted, however, that in this system the indirect 
3 Hc- 3 Hc interaction mediated by ripplons in the underly- 
ing 4 Hc film, which is not considered in most theoretical 
works, might be important. In addition, Nuclepore is be- 
lieved to be a much less uniform substrate than graphite. 

Recently, Sato et al. [ll[ found the G-L transition with 
p C Q ~ 1 nm~ 2 in the heat-capacity measurements on the 
third layer of 3 He on graphite down to T — 1 mK. This 
was inferred from a linear p— dependence of 7, the co- 
efficient of the leading T-linear term of C in the de- 
generate region, as well as a kink at 7 m Tideal- Here 
7ideai(= 7rfcfjAm/3?i 2 ) is the 7 value of an ideal Fermi 
gas spreading over the whole surface area (A) of the sub- 
strate, and m is the bare mass of 3 He. Note that 7 de- 



pends only on A and m not on the number of particles 
in the 2D case. One possible explanation for their result, 
which contradicts existing theory, is that, in the third 
layer, the relatively large plane-normal motion may sta- 
bilize the liquid phase (the quasi two-dimensionality). A 
variational calculation [12J supports this scenario but the 
subsequent ones [H, 0] do not. This hypothesis can be 
tested by extending their C measurement to the first or 
second layer in which the substrate confinement poten- 
tial is much deeper. The other issue is the role of surface 
heterogeneities in Grafoil [l3| , an exfoliated graphite sub- 
strate, used in most of the previous experiments includ- 
ing Ref. [llj]. This substrate is known to have a platelet 
(micro-crystallite) structure with a mosaic angle spread 
of about 30 degrees [13] and a platelet size of 10^100 nm 
15[. The role can be checked, for instance, by comparing 



results on the first layer of 3 Hc, which is directly on the 
Grafoil, and those in the upper layers. 

In this Letter, we report a result of new heat-capacity 
measurements of three different 3 Hc monolayers, i.e., the 
first, second, and third layers of 3 He on graphite, at very 
low densities never explored before, using the same ex- 
perimental setup as in Ref. [ll[ . We could determine the 
substrate heterogeneity effect explicitly in the first-layer 
measurement. By preplating the first layer with non- 
magnetic 4 He, the effect can thoroughly be removed in 
the second-layer measurement. Surprisingly, all the three 
monolayers show the G-L transitions with approximately 
the same p c o values (0.6 ~ 0.9 nm~ 2 ). This indicates that 
the quantum gas phase is not the ground sate of 3 He in 
strictly 2D and gives rise to a challenge for current many- 
body theories. 

In Fig. Ufa), we show measured C data for the first 
layer of 3 He adsorbed directly on the Grafoil substrate 
with A = 556 m 2 . This surface area is determined from 
the sub-step structure in N2 adsorption isotherm mea- 
surement corresponding to the x y3 commensurate 
phase formation in the first layer. Two heat-capacity 
contributions with distinct T-dependences develop suc- 
cessively as a function of p. The data at any densities 
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FIG. 1. (a) Heat capacities (C) of the first layer 3 He on 
Grafoil. The numbers are densities in nm" 2 , and those not 
denoted are 1.05, 1.25, 1.45 and 1.70 nm -2 , respectively, from 
bottom to top. The solid lines are fittings to Eq. ((TJ. The 
dashed line is C — 7idcair for a degenerate Fermi gas spread- 
ing over the whole surface, (b) Density dependence of the 
fitting parameter /3 in Eq. |T]). (c) Density dependence of 
the fitting parameter 7 in Eq. fTJ. The open circles are from 
Ref. [l(|. Only after the growth of /3 is saturated, 7 starts 
to increase linearly with p above 0.6 nm -2 . The horizontal 
dashed line represents 7 = 7idcai- 



above 0.3 nm 2 can be well fitted to 
C(T,p)^ 1 T-aT 2 + /3C a 



(1) 



in the T-rangc from 4 to 80 mK. The first two terms 
on the right-hand side of Eq. (JTJ are characteristic of a 
degenerate 2D Fermi liquid with spin fluctuations 17 1. 
Camor(T) is the heat capacity of the 0.45 nm" 2 sample. 
This is associated with nuclear spin degrees of freedom of 
amorphous 3 He 18 1 trapped on strong adsorption sites of 
Grafoil. The unusually weak T-dependence is a result of 
a wide distribution of exchange interaction in the amor- 
phous state. As shown in Fig. DJb), with increasing p, 
only the amorphous component (fitted /3 value) increases 
linearly with a small offset of 0.1 nm -2 , and is saturated 
above p « 0.6 nm" 2 . After then, the fitted 7 value starts 
to increase linearly until a kink at p = 1.4 nm -2 near 7 = 
7idcai- The density variation of 7 above 0.6 nm" 2 is very 
similar to that observed in the third layer above the in- 
tervening region (llj . To our knowledge, the only reason- 
able explanation for this is the phase separation between 
a degenerate Fermi liquid {puddles) with an almost fixed 
density of p c o (= 0.8 nm" 2 ) and a dilute gas phase with 
a negligibly small C contribution (see later discussion). 
It is clear that the G-L transition in the first layer devel- 
ops on the uniform region of the substrate independently 
of the preceding occupation of the heterogeneous sites by 
3 Hc. Such sites would be located only near platelet edges. 
The number of 3 He atoms contributing to C amor is 10% 
of that on the uniform surface to complete the v^xy^ 



commensurate phase [lj|. This ratio is consistent with 
previous thermodynamic measurements [20j . 

Next, we made heat capacity measurements of the sec- 
ond layer of 3 Hc on Grafoil prcplatcd with a monolayer of 
4 He, which preferentially occupies the first layer because 
of its smaller zero-point energy than 3 He. This tech- 
nique has widely been employed in previous experiments 
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2lL [22| . We introduced exactly the same amount of 
4 He (12.09 nm" 2 ) as that in Ref. [il|. As seen in Fig. 
[2ja), any non- Fermi liquid C contributions are absent 
here, indicating thorough elimination of the substrate 
heterogeneity effect by the 4 He prcplating. The data can 
be fitted to the formula: 



C(T) = 7T - qT 2 



(2) 



very well. The fitted 7 follows perfectly the p-linear de- 
pendence with a negligibly small offset (0.02 nm" 2 ) as 
well as a kink at p C Q = 0.6 nm" 2 and j c q = 1 ■S'Jideal 
(see Fig. [2{b)). Therefore, a G-L transition is observed 
again. Moreover, the second layer of 3 He should be the 
best representative of monolayer 3 He on graphite without 
heterogeneities. According to the previous experimental 
24 1 and theoretical Q determinations of the second layer 



promotion density of 4 He (11.4 ~ 11.8 nm -2 ), we expect 
that a small fraction (0.3 ~ 0.7 nm" 2 ) of 4 Hc is promoted 
to the second layer and preferentially occupies deeper po- 
tential sites above the substrate heterogeneities. This ex- 
plains why we don't observe C am or nor a sizable interven- 
ing region prior to the puddle region in the second-layer 
measurement. Our 7 data follow smoothly the previ- 
ous data [2^] using exactly the same experimental setup 
(open circles in Fig. E^b)) at p > p c0 , where 7, and hence 
the quasiparticle effective mass, increases progressively 
due to particle correlations. 

Let us briefly comment on possible finite-size effects 
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FIG. 2. (a) Heat capacities of the second layer of 3 He on 
Grafoil preplated with a monolayer of 4 He. The solid lines 
are fittings to Eq. ([2]). (b) Density variation of the fitting 
parameter 7 in Eq. (J2J. The open circles are from Matsumoto 
et al. [25|. Note that their data point at 2.38 nm -2 is not 
shown in Ref. [2f|. Otherwise, notation here is the same as 
in Fig. CD 
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caused by the platelet structure of Grafoil. The en- 
ergy discreteness estimated from the platelet size is 2 ~ 
200 /iK. This will not affect at least the leading 7T-terms 
in Eqs. ^ and and hence our G-L transition sce- 
nario, within the temperature range we studied (T > 
2 mK). On the other hand, the correction term ctT 2 , 



which is due to the spin fluctuations [17j], is suppressed 
depending on the puddle size within the two-phase coex- 
istence region of the second layer. Eventually, a/7 de- 
creases from 5 K _1 to zero with decreasing p presumably 
because of the long-wavelength cutoff of the fluctuations. 
More details of the size effects will be discussed elsewhere 



(a) 

100 

CM 



1 ; 1 1 
Puddle 


' 1 1 
C2+IC2 IC2 " 


Ilia 1Mb 






/ - 










*tr\ , 1 


third layer 



(b) 



b b fc> 6 



(c) 



(d) 



8 9 10 

P (1m" 2 ) 



FIG. 3. (a) Density variation of 7 deduced from measured 
heat capacities of the third-layer liquid 3 He on graphite pre- 
plated with a monolayer 4 He; present work (closed circles), 
Ref. [ll| (open circles). The magnetic contribution from the 
second-layer solid 3 He has already been subtracted as de- 
scribed in Ref. The solid line is guide for eye. The 
uniform liquid region above 8.1 nm -2 is divided into two re- 
gions, (C2 + IC2) and IC2, depending on the structure of 
the second-layer solid 3 He (see Ref. [jj]|). Schematic cross- 
sectional views of the first (b), second (c) and third- layer (d) 
3 He puddles on graphite. Only topmost three graphene-layers 
are drawn here. 

We have made additional heat-capacity measurements 
for the third layer of 3 He to understand further details of 
the density variation of 7 studied in Ref. [IlJ . Within the 
intervening region between 6.8 and 7.3 nm" 2 (Region Ilia 
in Fig. |3ta)), it was found that 7 varies in proportion to 
p with a factor of three smaller slope than that in the fol- 
lowing main puddle region of 7.3 < p < 8.1 nm -2 (Region 
Illb). We speculate that, in Region Ilia, promotion to 
the third layer as liquid puddles and compression of the 
second layer proceed simultaneously. This speculation 
is supported by the observed increases of magnetic C- 
isothcrms below 1 mK by about 10% in the corresponding 
density region [ill ]. A similar intervening re gion is also 



observed in the previous NMR experiment [22j, where 
the second layer of 3 He is compressed by adding 4 He. 
The compression of the second layer can either be so- 
lidification of a remnant high-density liquid, which may 
exist nearby the heterogeneities, or introduction of iter- 
stitial atoms to the commensurate phase (C2) [2^]. Con- 



sequently, we estimate p cu in the third layer as 0.9 nm" 2 
or slightly less. In Figs. [3fb), (c) and (d), cross-sectional 
views of the first, second and third-layer puddles of 3 He 
are imaged, respectively. 

The zinc superconducting heat-switch we used unfortu- 
nately does not allow us to extend our heat-capacity mea- 
surements beyond 80 mK where one expects to observe C 
anomalies associated with finite-T G-L transitions (T c ). 
We speculate that the highest T c (T c max ) is realized at 
P ~ Pco/2 m a T-range betwen 80 mK and 0.7 K. The 
high-T bound comes from the known T c max for 4 Hc in 2D 
Q . Then, a naive question is why we don't observe any C 
contributions from the phase-separated gas phase in the 
puddle regions. The dash-dotted line in Fig. 2] is a G- 
L phase separation line calculated for classical adatoms 
interacting with the Lennard- Jones potential [27| , where 
we adjusted the line so as to give /? cU = 0.6 nm -2 and 
j, max _ m j£ This T c max value was chosen arbitrar- 
ily. The low-density branch and high-density one give the 
equilibrium gas density p c a (T) and liquid one pj (T), 
respectively. Since they vary exponentially with T at 
T < T c max , we expect p c 9 /p c0 < 0.02 and p c l /p c0 « 1 
at T < 80 mK. That is why the C contribution from the 
gas phase is immeasurably small and the liquid phase is 
always degenerate with nearly the constant density p cu 
in our measurement. If the G-L phase separation does 
not occur, we should observe a smooth approach of 7 to 
7ideai when p decreases down to the lowest density sample 
(0.125 nm" 2 ; Fermi temperature (Tp) = 63 mK) without 
any kinks, which is of course totally not the experimental 
case. 
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FIG. 4. Low-density phase diagram of monolayer 3 He on 
graphite. The thick solid lines are the gas-liquid (G-L) transi- 
tion lines determined from this experiment. The dash-dotted 
line is a calculated one normalized to p c o = 0.6 nm" 2 and 

T raax = mK ( Ref _ daghed lme re p re sentS 

T F of 2D 3 He gas. The arrows denote sample densities at 
which we made the C measurements for the second layer of 
3 He. The inset shows a schematic top view of phase-separated 
liquid puddles in the second layer. 

The fact that the first three layers of 3 He on graphite 
have nearly the same p c o values (= 0.8, 0.6, 0.9 nm -2 ) ex- 
cludes the quasi two-dimensionality from possible expla- 
nations for the self-binding. This is because the confine- 
ment potentials and wave-function overlappings between 
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the successive layers are quite different each other in these 
monolayers. Furthermore, the indirect interaction (Vi n d) 
mediated by excitations in the underlayer should be quite 
different, too. Schick and Campbell [28| calculated V in d 
due to phonon exchange to be proportional to a factor 
n s c^ 2 e 2 in the case of substrates occupying a half-infinite 
space. Here n s , ct and e are the three-dimensional den- 
sity and phonon velocity of the substrate and the min- 
imum of the He-substrate potential, respectively. If we 
apply this theory to the present problem, this factor is, 
at least, an order of magnitude smaller for the first layer 
compared to the third layer. This is inconsistent with 
the fact that nearly the same p C Q values are obtained in 
these layers (2§|. We thus conclude that the observed 
G-L transition in the present experiment should be an 
intrinsic property of 3 He in strictly 2D. 

Our conclusion gives rise to a conflict with the ex- 
isting many-body calculations for 3 He in 2D 0-0 • It 
is, however, worthwhile to remark on recent variational 
and diffusion Monte Carlo calculations by Kilic and 
Vranjes [3(| on binding energies of 3 He molecules of N 
atoms in 2D. They obtained tiny but finite binding ener- 
gies, -(0.02-0.04) ml, for 2 <N< 6. Since the binding 
energy should decrease with N — > oo, their calculations 
seem to be consistent with the present experimental re- 
sult. 

In summary, we found the gas-liquid transition in the 
three different 3 He monolayer systems, i.e., the first, sec- 
ond and third layers, on graphite from the heat capacity 
measurements at low densities never explored before in 
the degenerate temperature region down to 2 mK. The 
phase-separated liquid phases have surprisingly similar 
densities (p c o = 0.6 ~ 0.9 nm -2 ) despite their quite dif- 
ferent environments, which indicates that 3 He atoms in a 
strictly 2D space arc self-bound forming liquid puddles at 
the ground state. The mean interatomic distance in this 
puddle is very large (1.1—1.4 nm). This would be to our 
knowledge the lowest density liquid ever found in nature. 
The present result contradicts the existing many-body 
calculations for 3 He in 2D, providing an important con- 
straint on theory. In future, it will be highly desirable to 
detect directly the expected thermodynamic anomalies 
at the critical temperature. 

We thank Tony Leggett for illuminating discussions 
and suggestions. This work was financially supported by 
Grant-in-Aid for Scientific Research on Priority Areas 
(No. 17071002) from MEXT, Japan and Scientific Re- 
search (A) (No. 22244042) from JSPS. DS acknowledges 
the support from JSPS Research program for Young Sci- 
entists. 



* |mroshi®phys.s.u-tokyo.a c.jp| 
[1] Note that laser-cooled alkali atoms stay gaseous down to 



the nanokelvin region as a metastable state only at very 
low densities where three-body collisions are not signifi- 
cant. 

[2] A. D. Novaco and C. E. Campbell, Phys. Rev. B 11, 2525 
(1975). 

[3] M. D. Miller and L. H. Nosanow, J. Low Temp. Phys. 

32, 145 (1978). 
[4] B. Krishnamachari and G. V. Chester, Phys. Rev. B 59, 

8852 (1999). 

[5] V. Grau, J. Boronat and J. Casulleras, Phys. Rev. Lett. 

89, 045301 (2002). 
[6] C. Urn, J. Kahng, Y. Kim, T. F. George, and L. N. 

Pandey, J. Low Temp. Phys. 107, 283 (1997). 
[7] D. S. Grey wall, Phys. Rev. B 41, 1842 (1990). 
[8] D. S. Grey wall, Phys. Rev. B 47, 309 (1993). 
[9] M. E. Pierce and E. Manousakis, Phys. Rev. B 62, 5228 

(2000). 

[10] B. K. Bhattacharyya and F. M. Gasparini, Phys. Rev. B 

31, 2719 (1985). 
[11] D. Sato, D. Tsuji, S. Takayoshi, K. Obata, T. Matsui, 

and H. Fukuyama, J. Low Temp. Phys. 158, 201 (2010). 
[12] B. Brami, F. Joly, and C. Lhuillier, J. Low Temp. Phys. 

94, 63 (1994). 

[13] Grafoil is a trademark of GrafTech International Ltd. 
[14] S. Takayoshi and H. Fukuyama, J. Low Temp. Phys. 158, 
672 (2010). 

[15] Y. Niimi, T. Matsui, H. Kambara, K. Tagami, M. 
Tsukada, and H. Fukuyama, Phys. Rev. B 73, 085421 
(2006). 

[16] D. S. Greywall and P. A. Busch, Phys. Rev. Lett. 65, 64 
(1990). 

[17] M. Ogura and H. Namaizawa, J. Phys. Soc. Jpn. 66, 3706 
(1997). 

[18] A. Golov and F. Pobell, Phys. Rev. B 53, 12647 (1996). 

[19] The 3 He density on the uniform region (p') can be dif- 
ferent from the nominal density (p) calibrated only at a 
single density, e.g., the density at the \/3 x \/3 com- 
mensurate phase in the first layer (6.37 nm -2 ). For 
our Grafoil substrate, the relation between p' and p is 
p'(nm~ 2 ) = 1.10(/9 — 0.6) at least for the 1st layer. 

[20] R. L. Elgin, J. M. Greif, and D. L. Goodstein, Phys. 
Rev. Lett. 41, 1723 (1978); K. Ishida, M. Morishita, K. 
Yawata, and H. Fukuyama, Phys. Rev. Lett. 79, 3451 
(1997). 

[21] C. P. Lusher, B. P. Cowan, and J. Saunders, Phys. Rev. 

Lett. 67, 2497 (1991). 
[22] E. Collin, S. Triqueneaux, R. Harakaly, M. Roger, C. 

Bauerle, Yu. M. Bunkov, and H. Godfrin, Phys. Rev. 

Lett. 86, 2447 (2001); E. Collin, Yu. M. Bunkov, and H. 

Godfrin, J. Phys.: Condens. Matter 16, S691 (2004). 
[23] T. Takagi, private communication. 

[24] M. Bretz, J. G. Dash, D. C. Hickernell, E. O. McLean, 
and O. E. Vilches, Phys. Rev. A 8, 1589 (1973). 

[25] Y. Matsumoto, D. Tsuji, S. Murakawa, H. Akisato, H. 
Kambara, and H. Fukuyama, J. Low Temp. Phys. 138, 
271 (2005). 

[26] D. Sato, K. Naruse, T. Matsui, and H. Fukuyama, sub- 
mitted. 

[27] S. Ostlund and A. N. Berker, Phys. Rev. B 21, 5410 
(1980). 

[28] M. Schick and C. E. Campbell, Phys. Rev. A 2, 1591 
(1970). 

[29] If we extrapolate the smooth density variations of 7 above 
the kinks to zero density, we have (1.2 ~ 1.3)7idcai for 



5 



the three monolayers. This is indicative of mass enhance- 
ment of a single 3 He atom caused by interactions with 
the underlayers. If so, Vind may play certain roles in the 
observed self-binding. Further theoretical studies in this 



aspect are desirable. 
[30] S. Kilic and L. Vranjes, J. Low Temp. Phys. 134, 713 
(2004). 



